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Despite the fact that the fundamental physical laws are symmetric in time, most observed pro- 
cesses do not show this symmetry. Especially the phenomenon of decay seems to involve a kind of 
irreversibility that makes the definition of a microscopic arrow of time possible. Such an intrinsic 
irreversibility is incorporated within the Rigged Hilbert Space quantum mechanics of the Brussels 
School, contrasting to the statements of standard quantum mechanics. As shown in this paper, the 
formalism bears significant advantages in the description of decaying systems, however the breaking 
of time symmetry can be avoided. 



I. INTRODUCTION 

Asking for the time evolution of a system means solv- 
ing the Schrodinger equation ifi^ip = Hip. Restricted to 
the case of closed systems - which is done from now on - 
this equation can formally be solved by 



(1) 



with ?i = 1 for the sake of simplicity and ip = ip{0)- The 
whole dynamics is thus governed by the unitary propaga- 
tor U{t) = e-'"* which maps any given "initial state" ip 
at time t = to its propagated counterpart tp{t) at time 
t. Quantum mechanical determinism is expressed by the 
fact that ip can also be regarded as a "final state" into 
which a state tp{—t) evolves by satisfying U{t)ilj{—t) — tp 
(see fig. 0). 
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FIG. 1. In standard quantum mechanics every state of a 
closed system can at the same time be regarded as an initial 
state and a final state. 

Since equation ^ is valid for all times i G R, giving 
tp{—t) = U{—t)ip, one is led to the operator identity 



U{~t)U{t) = 1 



(2) 



In other words: The set of evolution operators {U{i)} 
forms a unitary group U with the group property (g), 
giving any propagator U{t) its inverse element U{-~t) = 



W{t). The physical meaning of this mathematical state- 
ment is profound: in quantum mechanics there is no mi- 
croscopic arrow of time. Any state tp €H can be propa- 
gated to the future as well as to the past, so that there 
is a complete time symmetry. 

In order to explain an arrow of time despite these facts, 
there are several ideas. A few of them are: 

1. Any measurement on the system interrupts the de- 
terministic evolution of the state by applying ran- 
dom projections (plus normalization) instead of 
unitary operations (the "collapse of the wave func- 
tion"). This process cannot be reversed. 

2. Real physical systems can never be regarded as 
closed. There is always an interacting environment 
with many unobservable degrees of freedom, induc- 
ing non-unitary time evolution of the observed sys- 
tem. 

3. Time symmetry in classical mechanics is only bro- 
ken with respect to macroscopic properties of a 
system (e.g. temperature or entropy). Consider- 
ing macroscopic observables in quantum mechan- 
ics, this effect also shows up (coarse graining argu- 
ment). 

All these ideas have been worked out to elaborate phys- 
ical models with strong experimental support and all of 
them break time symmetry (see e.g. [[l5|j9|). So one could 
be content with the situation and regard those models as 
fairly good explanations of an arrow of time, which is 
indeed an undeniable part of human experience. 

A group of scientists around Ilja Prigogine (nobel prize 
in chemistry 1977), referred to as the "Brussels School", 
is not satisfied with the situation and has worked out 
a quantum mechanical formalism especially applying to 
unstable systems (i.e. systems involving decay), where 
time symmetry is claimed to be broken at the micro- 
scopic level (see e.g. ll|,|,|||,|,0). 
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FIG. 2. In Brussels School's quantum mechanics every 
state is whether an initial or a final state. Time evolution 
falls apart into two different directions. 

This fact is expressed by the splitting of the unitary 
group U into two semigroups U- and L(+ of propagators 
for negative and positive times, such that past and future 
states can no longer be regarded as equivalent counter- 
parts. This leads to a so-called "intrinsic irreversibility", 
implying a microscopic arrow of time (see fig. 

This paper is separated into two different parts. The 
first part is dedicated to a review of the Rigged Hilbert 
Space formalism worked out by the Brussels School. It is 
shown how it works, how it applies to physical systems 
and what its advantages over standard quantum mechan- 
ics are. The second part is devoted to the role of time 
symmetry within the Rigged Hilbert Space formalism. 
Apart from its practical advantages, the Brussels School 
points out that a microscopic arrow of time can be im- 
plemented here. It will be shown that such an implemen- 
tation is possible but not mathematically necessary. 



where H\E,i) = En\En), H\E) = E\E) and p{E) is the 
energy density. As soon as the energy decomposition of 
the initial state generally given by 

IV-) =^V^„|i?„)+ / dEp{E)i,{E)\E) , (4) 

is known, the propagated state is also known and reads 
|^(i))=^V'.e-'^"*|i?„) 

n 

+ I dEp{E)ij{E)e-''''\E) . (5) 

The expectation value of an observable A can be calcu- 
lated for any time t using 

<A>{t) = m)\Am)) ■ (6) 

So the fundamental key to the dynamics of a system is 
the spectral decomposition of its Hamiltonian. All dy- 
namical properties are monitored by this decomposition 
and here is where the new formalism enters. There is 
a class of systems where the Brussels Schools's spectral 
decomposition of the Hamiltonian differs from the usual 
decomposition: systems envolving decay, called unstable 
systems. 



B. Resolvent techniques 



II. THE BASIC CONCEPTS 

A. Spectral decompositions 

Writing down the formal solution (|^) of a closed sys- 
tem's Schrodinger equation is not yet really satisfying. 
How can expectation values of observables be numeri- 
cally calculated? To this aim one has to find a spectral 
decomposition of the system's Hamiltonian iJ, thus one 
has to solve the time-independent Schrodinger equation 
Hil> — Eip. Since iJ is a self-adjoint operator, its spec- 
tral decomposition is complete and orthogonal and may 
therefore be used as a spectral decomposition of the unity 
operator 1. Furthermore the spectrum of H is real and 
consits in general of a discrete part and a continu- 
ous part ac- In special cases there are other parts, e.g. 
singular continuous and dense spectra, but these cases 
are excluded here. With this restriction, using the Dirac 
bra-ket notation, the spectral decomposition of the unity 
operator reads in general 

l = Y,\En){En\+ f dEp{E)\E){E\ , (3) 



The basic concept of the Brussels School's formalism is 
to find a different spectral decomposition of the Hamil- 
tonian. In contrast to the method of complex scaling 
(see Q), the Hamiltonian is not modified and remains 
self-adjoint; it only acquires a different representation in- 
volving generalized eigenvectors contained in an extended 
distribution space. The modified decomposition can most 
easily be achieved by using the Hamiltonian's resolvent 



The spectrum of H is defined by the irregularities of its 
resolvent. The discrete eigenvalues are the singularities 
of R{z), i.e. those points in the complex plane where 
the operator R is not everywhere defined on the Hilbert 
space Ti. The continuous eigenvalues are those points 
where R is an unbound (and thus non-continuous) op- 
erator on Ti. Anywhere else on the complex plane the 
resolvent is an analytic function mapping complex num- 
bers onto the set of linear bounded operators. As H is 
self-adjoint, any irregularities of R take place on the real 
axis, where the discrete eigenvalues form the point spec- 
trum and the continuous eigenvalues define a cut in the 
complex plane (see fig. §). 
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FIG. 3. The eigenvalues of the Hamihonian H are the irreg- 
ularities of its resolvent R{z) = . Since H is self-adjoint, 
they all take place on the real axis, where the discrete eigen- 
values En G (Jd are singularities and the continuous eigenval- 
ues E £ ac form a cut in the complex plane. 

It is possible to find an analytic continuation of R from 
both sides of the cut, each leading to a different operator 



R±{E) 



1 



H±ie 



rUe), 



(8) 



with e +0 understood here and in the following. Note 
that this continuation is not possible at singularities em- 
bedded in the continuous spectrum. Note further that on 
the real axis outside the continuous spectrum both op- 
erators i?+ and i?_ coincide with R. The two continued 
resolvents are known as the retarded and the advanced 
Green operator, respectively. 

The resolvent can be used to construct eigenvectors. 
The residuum of the resolvent at a discrete eigenvalue 
yields the dyadic product of the corresponding eigenvec- 
tor: 



lim (z 



EMz) = \Et){Ei 



(9) 



Note that both dyadic products project onto the same 
subspace corresponding to En- If the discrete eigenvalue 
is isolated, i.e. not embedded in the continuous spec- 
trum, both limits from above and below the real axis co- 
incide. In the following we will assume that all discrete 
eigenvalues are isolated and we therefore have 



lim {z~En)R{z) - \E„){En\ 



(10) 



A discrete eigenvector |£'„) of H can be obtained by ap- 
plying the residual resolvent to any Hilbert vector lip) 
being not orthogonal to the eigenspace of : 



\Er,) = lim {z - E„)R{z)\i:) 



(11) 



where ip„ — (Enlip) ^ is an adequate normalization 
constant, such that {En\Em) — Snm- 

In analogy to the discrete case the residuum of the re- 
solvent at one of the continuous eigenvalues yields the 



dyadic product of the corresponding continuous eigen- 
vector: 



Mm {z-E)R{z) = \E^){E^\ 

z—tE±ie 



(12) 



Note that both dyadic products are equivalent represen- 
tations of the 5-Operator 6{E — He), defined by 

/ dEp{E)S{E-H,)f{E):= f{H,) , (13) 



where He is the continuous part of the total Hamilto- 
nian H. Since the continuous eigenvalues are not iso- 
lated, there are always pairs of eigenvectors. They cor- 
respond to retarded and advanced solutions of the con- 
tinuous eigenvalue equation and are known as incoming 
and outgoing waves. Both solutions can equivalently be 
used in the spectral decomposition and they can be con- 
structed by 



\E^) 



1 ±ie 



(14) 



tP±{E) E - H ±ie' 

where Tp±{E) = {E^\ip) ^ is a suitable normalization 
constant, such that {E^\E'^) = 5{E - E'). 

Using dyadic products, the retarded and advanced 
Green operators can be written as 



Eea,: R±{E) 



Ve 

E-H 



TiAE){E\ , (15) 



where \E){E\ is any of the two dyadic products \E^){E'^ \ 
and the principal value distribution Ve is defined by 



r pE+e 

/ dE'VEf{E') := / dE' f{E 



(16) 



The difference of the two Green operators gives the 
dyadic product of the continuous eigenvector: 



-{R+{E)-R_{E)}^\E^){E^\ 



(17) 



Altogether, the spectral decomposition of the unity op- 
erator may be expressed by the exclusive use of the re- 
solvent: 

1 = ^ lim (z - En)R{z) 

dE p{E) {R+{E) ^ R^{E)] . (18) 



'2^ 



This unity decomposition will be used as the starting 
point of the Brussels Schools's unity spectral decomposi- 
tion. 
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C. The second Riemann sheet 

The resolvent can be analytically continued to the com- 
plex plane beyond the real axis entering the second Rie- 
mann sheet: 



R±{z) lim R±{E), 



E eac 



(19) 



On the one half of the complex plane where the contin- 
uation started (the first Riemann sheet), both operators 
coincide with the original resolvent. On the other half 
(the second Riemann sheet) it might happen that new 
singularities appear, and that in a pairwise manner: If 
Zk is a singularity of R+{z) in the lower half of the com- 
plex plane, then z^ is a singularity of R^{z) in the upper 
half (see fig. U). Unstable systems are exactly those with 
a resolvent having complex poles on the second Riemann 
sheet. 

RM) 

2nd Riemann sheet 




Re^ 



FIG. 4. The second Riemann sheet. New singularities of 
the analytically continued resolvent may appear, both com- 
plex conjugated to another. 

The second Riemann sheet poles act as complex reso- 
nances and cause the decay of certain states. 



D. An extended spectral decomposition 

Consider the continuous part of the unitary decompo- 
sition (|l|): 



2^ 



/ dE p{E) {R+{E) - R^{E)] 

J (7 c 



(20) 



The following manipulations of Ic are reducing its do- 
main. As a consequence one is forced to use a smaller test 
space, enlarging at the same time the distribution space 
and inducing this way a certain Gelfand tripel or Rigged 
Hilhert Space where the modified completeness relation 
remains valid. We delay this discussion to section 

If one analytically deforms the integration path along 
the continuous spectrum Oc to a path r+ leading through 
the lower half of the complex plane, one enters the second 



Riemann sheet and crosses the singularities of i?+ (z) . Let 
the set of these singularities be defined by cr+ = {zk\ and 
be enclosed by ctc U r+. In the course of the deformation 
circular integrations around the singularities are split off, 
each leading to a residuum of R+{z) at these points (see 
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FIG. 5. Analytically deforming the integration path along 
Oc to a path r+ one crosses the singularities of R+(z), which 
have to be split off from the integration. 

The modified operator now reads 

1+=^ hm {E-Zk)R+{E) 

k 

+ 1. f dzp{z){R+{z)~R^iz)} . (21) 
Introducing a complex distribution D±{E), defined by 
dED±{E)f{E) := / dzf{z) , (22) 

where r_ = T*^ and /(z) being the analytic continua- 
tion of / to respectively the upper and lower half of the 
complex plane, one can rewrite as 

1+=^ hm {E~Zk)R+{E) 

k 



2tt 



dE p{E) D+{E){R+{E) - R^{E)] .(23) 



In analogy to the standard case one defines dyadic prod- 
ucts of generalized eigenvectors of H by 



!/+)(/,- 1 hm {E-Zk)R+{E) 

for the discrete case and, using (|l7|), 

\f+{E)){r{E)\=D+{E)\E){E\ 



(24) 



(25) 



for the continuous case, where \E){E\ is any of the two 
dyadic products \E^){E^\. Hence we have 
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/ dEpiE)\E){E\=J2\f^){fk 



dEp{E)\f+{E)){r{E)\. (26) 



The generalized eigenvectors appear as pairs of right and 
left eigenvectors. The discrete eigenvectors are called 
Gamov vectors and can be constructed by: 



1/+)- — lim {E~Zk)R+{Em, 

iPk E^Zk 



(.4^1 = — 1™ {E~z,;){m+{E), 

Vk E^Zk 



(27) 
(28) 



where is a suitable vector and V-'fc = {fk\i^) 7^ is 
a normalization constant, such that {f^\f^) = Ski- By 
adjunction one sees that \f^) ^ IZ/T)- Furthermore we 
have if^lf^) = WfkW " ^- These objects are left and 
right eigenvectors of the Haniiltonian H corresponding 
to complex eigenvalues: 



{fk\H^{fkW 



(29) 
(30) 



Since the Gamov vectors have zero norm, they all co- 
incide in Hilbert space with the null vector, which is 
not allowed to be an eigenvector. Also, since H is self- 
adjoint, it cannot have any eigenvectors corresponding 
to complex eigenvalues. Concluding, the Gamov vectors 
cannot be elements of the Hilbert space. In fact they 
are distributions in a suitably chosen distribution space 
(see section IVB). The continuous eigenvectors can be 
constructed in two ways: 

|/±(i?)) = D^{E)\E+) ^ {f^{E)\ = D^{E+l (31) 
\g^{E)) ^ D±{E)\E-) ^ {g^{E)\ ^ D^{E){E-\. (32) 

By application to test vectors with analytic extensions to 
the complex plane one sees that 

D+{E)\E-) = \E-) (33) 
D^{E)\E+) = \E+) . (34) 

Hence there are two equivalent bases given by 

B+ = {\Er,){E,-,l\f+){f^\,\f+{E)){E+\}, (35) 
B+ = {\E„){E„\, \E-){g-{E)\}. (36) 

Both bases can equally be used to span the Hilbert space 
Ti,. Arbitrarily we select and obtain a retarded unity 
decomposition 



U=J2\E,,){E,,\+J2\fk){f, 



k I' 



dEpiE)\f+iE)){E+\ 



(37) 



In complete analogy one can modify the spectral de- 
composition by deformation of the integration path to 
r_ = r'j_ in the upper half of the complex plane, enter- 
ing the domain of R^{z) with singularities at cr_ — {z^}. 
One then obtains the two bases 



such that 

\fk) 



{ft' 



B2 = {Bt)\ 



— Yym^{E - zl)R^{Em 

tpk E^zl 

^ liTn^iE - zimR^{E) 



(38) 
(39) 



(40) 
(41) 



^' ' i'k E- 

iriE)) ^ \E+), {f+{E)\ DAE){E+l (42) 



\g-{E))^D-{E)\E-), {g+{E)\^{E-\ 



(43) 



Arbitrarily choosing B^ one obtains an advanced unity 
decomposition: 



i_ = 5]|i?„)(£;„| + ^|/,7)(/^ 



dEp{E)\E+){f+{E)\ 



(44) 



which is just the adjoint of the retarded unity decom- 
position, so (1+)^ = 1_. Note, however, that both de- 
compositions are equivalent representations of the unity 
operator, thus in the sense of an operator identity we 
have 1 , = 1_ = 1. 



E. Using the extended decomposition 

The unity decompositions obtained in the preceding 
section can be used for different representations of vec- 
tors and operators. Using the retarded decomposition 
the Hamiltonian reads 



dEp{E)E\f+{E)){E+\ 



(45) 



The adjoint of the above expression yields the Hamilto- 
nian in the advanced decomposition. Since both decom- 
positions are equivalent representations, one sees that 
= 77, i.e. the Hamiltonian remains self-adjoint. 
The propagator in the retarded decomposition reads 

c/(i) = E^"^''"*l^")(^"l + E^"'''*IAt)(/fc"l 

n k 

+ f dEp{E)e-^'''\f+{E)){E+\ . (46) 

Since the complex numbers Zk are located in the lower 
half of the complex plane, the contribution from the 
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Gamov vectors decays in the future. If an initial Hilbert 
state has Gamov components, they will disappear for 
t —* oo. So if the state has no rotating components, i.e. 

— 0, the whole state will disappear in the future 
with the decay rates of the Gamov contributions given 
by the imaginary part of the contributing z^: 







(47) 



with il}^ = (/^IV'), Zk = Vk - ilk and Vk.ik > 0. Each 
contribution from a Gamov state also induces a certain 
energy shift, given by the real part of the Zk- The advan- 
tage of the extended spectral decomposition is getting 
clear now: Not only the stable frequencies can imme- 
diately be read from the spectrum, but the decay rates 
and the energy shifts, too. By additionally checking the 
amplitudes between the initial state and the Gamov 
vectors one has immediate access to the dissipative dy- 
namics of the system. 



F. A closer look 

Apart from the contribution of the Gamov vectors one 
has to consider the contribution from the generalized con- 
tinuous eigenvectors. Let us call it the background con- 
tribution. Since the decomposition involves the complex 
distribution \ f^{E)) = D+{E)\E) , a test vector (p has to 
be introduced. The background contribution then reads: 



(^l^bg(t)) :- / dEp{E)e-^^'{ip\f+iE)){E+\^) (48) 

= / dEp{E)D+{E)i;+iE)e-'^*^X{E) (49) 

where ?/'+(£') = {E+\ip) and (pX{E) = {'p\E+). Due to 
{E) the integration path is transformed into the curve 
r+ below the complex eigenvalues Zk- 



dzp{z)e'"'-''ijj{z)ip*{z) 



(50) 



where tp+{z) and ^p*i^{z) are the analytic continuations of 
respectively ip+{E) and (p'^{E) to the lower half of the 
complex plane (provided their existence). The curve r_|_ 
is arbitrary as long as r+ U CTc surrounds cr+, i.e., r+ is 
placed below the poles Zk, which all have negative imag- 
inary parts. For t > 0, due to the exponential factor 
e~*^*, the background contribution thus becomes more 
and more neglectable in comparison with the pole con- 
tributions. For initial states with high decay rates, i.e. 
pole contributions with small imaginary parts, the back- 
ground contribution may be neglected already for small 
times t > 0. This approximation is called the pole approx- 
imation and it is essentially equivalent to the Wigner- 
Weisskopf approximation (see [|l4| ). The time evolution 
of the initial state may be approximated by an exponen- 
tial decay and its energy spectrum by a Breit-Wigner 



distribution. This approximation is very popular, e.g. in 
quantum optics, and indeed in most cases it is a very 
good one. 



G. Do divergences break time symmetry? 

Obviously the above argument does not hold for neg- 
ative times. Here, the contributions from the Gamov 
vectors diverge. In several publications of the Brussels 
School (see e.g. or A. Bohm in ||^) this is an ar- 
gument in favour of a broken time symmetry, since for 
negative times the advanced decomposition ( p[ ) is the 
better choice. However, as will be shown below, this is a 
rather intuitive reasoning and does not enforce a break 
of time symmetry. 

Let us see what happens with the retarded decompo- 
sition of a Hilbert state for negative times. Of course the 
Gamov contributions now diverge, but they are compen- 
sated by the background contribution. For t —oo the 
background integral (|^ gains influence, since the inte- 
gration path r+ is placed below the complex eigenvalues 
Zk- The same argument, that for increasing t > the 
background contribution can be neglected in comparison 
with the pole contributions, is turned upside down for 
t < 0. The smaller imaginary part of r+ induces a higher 
divergence than the contribution of the poles. By con- 
struction their divergences cancel each other, since the 
sum of both yields the contribution of the standard con- 
tinuous eigenvectors (see eq. (|2^)), so a divergence of an 
initial Hilbert state backwards in time is prevented. The 
divergence of the Gamov contributions is even more than 
compensated: If a Hilbert state decays in the future, it 
will also decay in the past. To see this one can use the 



advanced decomposition (4^) to investigate the system's 
evolution to the past and one would find a completely 
symmetric situation with majorizing pole contributions 
decaying backwards in time. 

The situation is different if a Gamov vector \f^) is 
taken as the initial state. Now there is only one pole con- 
tribution Zk decaying in the future and diverging in the 
past. But Gamov vectors are only "existing" in the spec- 
tral representation of Hilbert states and they themselves 
do not refer to elements of physical reality, similiar to 
all other "generalized eigenstates" appearing in spectral 
decompositions. As constituted by the axiomatic founda- 
tions of quantum mechanics a physically realizable state 
is represented by a Hilbert vector (actually a family of 
Hilbert vectors normalized to unity and differing only in 
phase). Gamov states are no Hilbert vectors, so their 
properties are not expected to be physically meaningful. 

Concluding, both decompositions are equivalent, but 
the use of the retarded decomposition is more intuitive 
for t > 0. However, this should not be taken as the 
reason for a broken time symmetry. A more convincing 
argument in favour of a broken time symmetry is the 
splitting of test spaces. As we have seen in (bw , analytic 
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continuations of the test functions are needed to have 
the decompositions well-defined. In section IV we will 
get back to this point. 



\Un) = \UJ„) + 



A 



En — Ho 
LOn + \{uJn\W\Un) 



}7iW\Un) 



(56) 
(57) 



H. Perturbation theory 



Besides the advantages pointed out in section HE we 



now face another important advantage of the Brussels 
School's formalism. Perturbation theory can be applied 
even in cases where standard quantum mechanics fails. 
Those cases are featured by the phenomenon of insta- 
bility, which is strongly connected to the sensitivity of 
a system towards small perturbations. Let the system's 
Hamiltonian be of the form 



Hx = Ho + XW 



(51) 



where W is the perturbation operator, switched on by 
the perturbation parameter A G [0, 1], and Hq is the un- 
perturbed part of the Hamiltonian, whose spectral de- 
composition is known and reads 



Ho = ^UJn\^n){L^n\ + / dcj Pq{u;) Uj\u)) {uj 



(52) 



To simplify the discussion let the spectrum of Hq be non- 
degenerate, so LOn 7^ (^m for n ^ m. As long as the per- 
turbation strength A is not explicetely needed, it is set to 
1 and i?A=i is written as H. The spectral decomposition 
of H is yet unknown and reads in general 



H = J2En\Un){Un\ + / 



duj p{uj) Uli\u {uj)){u (a;)|. 



(53) 



where the sum is taken over the discrete spectrum dd and 
the integration is performed over the continuous spec- 
trum tJc oi H. 

The continuous spectrum Uc is not affected by the per- 
turbation (provided the perturbation operator W is Hq- 
compact, which is assumed here) and there is an identity 
mapping of perturbed and unperturbed continuous eigen- 
values, so in both cases the integration is performed over 
tJc- The discrete case is different. Switching to Brioullin- 
Wigner perturbation theory let us define pairs of orthog- 
onal projection operators 



|w„)(w„| and 
1-Pn ■ 



(54) 
(55) 



Furthermore let the perturbed eigenvector of H be 
not normalized to 1 but rather let (M„|ci;„) = 1, such 
that Pn\un) = Applying these projections to the 

perturbed eigenvalue equation H\un) one obtains the 
Brioullin-Wigner formulas for the shifted eigenstate |zi„) 
and the shifted energy i?„: 



Iteration leads to the corresponding perturbation series 
A 



En — Ho 



p=0 

oo f 

En = LUn + X^{uJn\W I Qn 

p=0 



QnW} \iUn), 



A 



En — Ho 



(58) 



nW\ \Un). (59) 



Take a closer look at the above formulas ( |58| ) and (|5^ 
They are only valid if the expression 



En — Ho 



ni^n 



En ~ l-^n 



+ / duj p{uj) 



En - LO 



(60) 



is well-defined. Here resonance appears as a major ob- 
stacle. There are two different types of resonance making 
a perturbation expansion of eigenstates and eigenvalues 
impossible: 

1. discrete resonance: The eigenvalue a;„ may be 
shifted in the course of the perturbation to a neigh- 
boured eigenvalue of Ho (see fig. ^). Expres- 
sion (|6^) is no longer defined and the corresponding 
perturbation series diverges on account of the sum- 
mation term. 

2. continuous resonance: If the eigenvalue w„ lies 
next to the continuous spectrum of Ho, it may be 
shifted there (see fig. ||b). Now the integral term 
in (|o|) causes a divergence of the perturbation se- 
ries. 

In both cases the perturbation expansion is limited to 
a finite convergence radius of the perturbation parame- 
ter A . Outside its range the shifted eigenvalue disap- 
pears from the spectrum of H\ and the mapping from 
the unperturbed to the perturbed spectrum is no longer 
one-to-one (see fig. |^) . 



b) 



En 



FIG. 6. a) discrete resonance, b) continuous resonance 
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Now regard the case of discrete unperturbed eigenval- 
ues embedded in the continuous spectrum. Here the con- 
vergence radius of A reduces to zero, i.e. the perturbation 
series are no longer analytic in A. An arbitralily small 
perturbation induces continuous resonance. As a conse- 
quence embedded states get unstable and start to decay. 
Decaying states cannot be part of the spectral decompo- 
sition of Hx, they disappear from there and cannot be 
reconstructed for A — > 0, i.e. Hx 7^ Hq. Hence pertur- 
bation theory is not aplicable. 



^^0 '-Hx 




or its iteration, the Born series: 

Take a closer look at the construction of Gamov states 
in eq. ( |6l] ) . The divergences due to continuous resonance 
are prevented by introducing a small complex shift in 
the denominator and analytic continuation to the second 
Riemann sheet. As a consequence, the energy eigenval- 
ues are shifted to respectively the lower and upper com- 
plex plane, inducing exponential decay and Brett- Wigner 
energy distribution of the corresponding Gamov states. 
The expressions are still analytic in A, so perturbation 
theory can be applied. This is a major advantage of the 
Brussels School's formalism. 



FIG. 7. In the course of the perturbation the continuous 
eigenvalues remain unaffected, while the discrete eigenvalues 
may disappear from the spectrum. In this case the mapping 
between the discrete spectra is no longer one-to-one and the 
perturbation series cannot be analytic. 

There is a way out: One has to include unstable states, 
i.e. Gamov vectors, into the spectral decomposition of 
Hx- In terms of perturbation theory the Gamov vectors 
can be constructed by modified BriouUin- Wigner equa- 
tions: 



|/f)-|c.,)+ lim Q,— 
z±=^,. + AK|M/|/±) 



A 



Ho±ie 



QkW\f^) 



(61) 
(62) 



where = are the singularities of i?+ on the lower 
half of the complex plane and z^T = z'^ are those of i?_ on 
the upper half. There is a one-to-one mapping between 
the disappearing discrete eigenvalues Uk and the com- 
plex eigenvalues Zk- Iteration leads to the corresponding 
perurbation series 



00 f 

|/±)= lim^J] Ofc- 



-^-—-QkwYliu,), (63) 

z^^u;k + X{u;k\W\f^) . (64) 
The generalized continuous eigenstates are 



\f+ico)) = D+{co)\u+iu)) 
{rii^)\ = {u+ico)\ , 



(65) 
(66) 



with the complex distribution defined in (|22|). The 
standard continous eigenvectors can be constructed using 
the Lippmann-Schwinger equation 



\u+iL0)) = \cu) + 



A 



Ho + ie 



W\u+iu;)) 



(67) 



III. EXAMPLE: THE FRIEDRICHS MODEL 

The Friedrichs model (K. Friedrichs 1948, see j|]) is 
the simplest model involving instability. Here one single 
discrete state is coupled to a continuum of states by an 
interaction Hamiltonian W. A physical analogon to this 
model is e.g. the Wigner- Weisskopf model (Wigner and 
Weisskopf 1930, see [Q) of a twolevel atom coupled to 
the vacuum field. Another example is the Auger effect of 
autoionisation (see e.g. ||ll|). 

The total Hamiltonian is of the form 



H = Ho + W 



(69) 



poo 

where Ho = uji\l){l\+ du uj \lu) {uj\ , (70) 
^0 

duji^W{uj)\uj){l\ + W*{uj)\l){uj\} . (71) 



Using the second resolvent identity 
R = Ro + RWRo 



(72) 



with the resolvent Rq(z) = —Kr- of the free Hamiltonian 
Ho and the resolvent R{z) 



z~H„ 

= of the total Hamilto- 

z — ri 

nian. The resolvent R{z) can be obtained by applying it 
to the free states |1} and \uj) and doing some rearrange- 
ment: 



R{z) 



1 



1 



du — ^ — ^jtij) 



duj ^ — - (uj 



Z — UJ 



duj 

z-uj 



uj){uj\ 



The complex function 

77(2) = Z - UJl 



duj 



\W{u)[^ 



(73) 



(74) 
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has analytic extensions 77±(w) = rjluj ± ie) from above 
and below to the real axis. The zeros of r/ on the second 
Riemann sheet define the complex poles of the resolvent: 



r;+(zi) = 0, 



(75) 



The pole zi can be obtained by numerically calculating 
the complex zero of 77+ (z). In first order one gets 



(76) 



The Gamov vectors can be obtained using eqns. (|2 
and (|2~ 



1/1+)-- lini {E ~ z,)R+m^) 

Ci E^zi 



duj 



[zi - w] + 



(/fH- hm {E - z^){l\R+{E) 

Ci E^zi 



auj — {uj 

[zi + 



where the complex distribution is defined by 



dhJ 



ip{uj) 



:= lim / duj ■ 

zi — wj+ £^21 Jq E — uj + ie 



(77) 
(78) 
(79) 
(80) 

(81) 



The standard outgoing states can be constructed using 
eq. (0): 



\u+{u;)) 



1 



c(cj) E- H + ie 
\uj) + ^ X 



|1) 



(82) 



(83) 



Using eq. ( pl|) the generalized continuous eigenvectors 
read 



\UiL0))=D+{L0)\u+iu)) 



(r(c^)i = (^+(c.)i 



Lo — Lo' + ie 



, , W{uj) 
(uj + — V4 X 



UJ — uj' — le 



(84) 

(85) 
(86) 

(87) 



where the complex distribution is defined by (|22 
with the curve r+ leading from to cx) below zi. 



In the pole approximation, i.e. neglecting the back- 
ground contribution from the continuous eigenvectors, 
the initial state |1) decays with a decay rate F given by 
two times the imaginary part of zi , 



p^{t) = \{i\um)\' 



-rt 



where T = 2 * Imzi « 27r|W^(wi)p, (89) 
which is in agreement with Fermi's Golden Rule. 

IV. TIME SYMMETRY 

A. The Rigged Hilbert Space 

When dealing with infinite dimensions, the Hilbert 
space appears to be too small in the following sense: Po- 
sition and momentum operator have no Hilbert eigenvec- 
tors and Dirac's bra-ket formalism is mathematically not 
justified. To correct this problem, the Hilbert space must 
be enlarged to include improper states like ^-functions 
and plane waves. The framework to enable these fea- 
tures is given by the theory of the Gelfand triplet, also 
called the Rigged Hilbert Space (RHS). 




FIG. 8. To deal with improper states like (5-functions and 
plane waves and to make Dirac's bra-ket formalism work, one 
has to leave the Hilbert space. 

The concept is the following (comp. e.g. p). Take 
the Hilbert space TC = oi squared integrable functions 
and find a suitable test space $ C 7i with the properties 

1. $ is a locally convex topological vector space, 

2. $ is complete with respect to its own topology, 

3. $ is dense in H, 

such that the topological dual of $, consisting of all 
linear continuous functionals on <I>, called distributions, 
includes the desired improper states. The Hilbert space 
itself is his own topological dual, i.e. Ti^ — H. For any 
test space $ being smaller than H, the Gelfand triplet 



$ C H C $ 



(90) 
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makes it possible to extend the bracket {x\f) to all pairs 
</? e X e by 



(91) 
(92) 



On pairs of Hilbert vectors the bracket conicides with the 
scalar product in Ti. The "kets" and "bras" are now iden- 
tified by linear functionals on the opposite dual space, i.e. 



and 



10 
(0 



(•10 
(0-) 



(93) 
(94) 



where ^ S W, $ or <i>^^. An operator A can be extended 
to the distribution space usme; 



(95) 



Hence A can only be extended to $t if its adjoint A^ is 
defined on and does not lead out of the test space <&: 



(96) 



The test space of standard quantum mechanics is the 
Schwartz space S of rapidly decreasing functions, defined 
by 



S:={f e C°°|sup|2;"9^7(2;)| < oo Vn,™} 



(97) 



The Schwartz space is a physically convenient test space, 
since here position and momentum operator are every- 
where defined, they fulfill (|9^), they are self-adjoint and 
have their eigenvectors among the tempered distributions 
in iSt. Furthermore, they keep their real spectrum, which 
makes the RHS (5, ^S^) be called a tight rigged Hilbert 
space. Another popular test space is the space V of 
C°°-functions with compact support. But V is not tight 
rigged, since the momentum operator obtains a complex 
spectrum with complex plane waves e*^^ G 2?^ as eigen- 
vectors. 



B. Breaking time symmetry 



Now we turn to the RHS of the Brussels School. As 



already stated in section IIP the definition of generalized 
eigenvectors is only possible if the test functions can be 
analytically continued to whether the lower or the up- 
per half of the complex plane. The retarded right eigen- 
vectors 1/^^) and |/+(i?)} involve complex distributions 
continuating the test functions <f*{E) = {f\E) to the 
lower half of the complex plane, denoted by (D-. There 
is a class of functions fulfilling this condition: the Hardy 
class from below (for detailed information on Hardy 
classes see 0). If (p*{E) is in then its complex conju- 
gate (p{E) = {E\if) is in the Hardy class from above, 
whose members can be analytically continuated to the 
upper half of the complex plane, denoted by (D+. In de- 
tail: The function (p{E) is in if and only if it has an 



analytic continuation to (D±, such that there is a C < oo 
and 



/oo 
dE\ip{E±iy)\^ < C 
-oo 



(98) 



Thus, restricted to the real line, Hardy functions are L^- 
functions. An equivalent definition can be given using 
the Paley-Wiener theorem: The Hardy spaces H"^ are 
formed by the inverse Fourier transforms of L^-functions 
with support on respectively the positive and negative 
semiaxis ]R±: 



Hi 



F-i{i2(R±)} 



(99) 



For physical convenience it is necessary, that the test 
functions are also in the Schwartz class S. Furthermore, 
their energy representation can be restricted to the pos- 
itive semiaxis, since energy is bounded from below by 
zero. Altogether, the test spaces of the Brussels School 
are defined by 



± •- 



{Hlns)\ 



The topological duals (E'j- then include the complex dis- 
tributions used in the generalized eigenvectors. The for- 
malism of the Brussels School can be applied in the RHS 



]R+ 



(100) 



(101) 



Since there are two RHS's, one has to investigate how to 
deal with them. The retarded unity decomposition ( |37| ) 
contains the dyadic products \fk){fk\- ^PPlyi^ig test 
vectors Lp+ and (f^ to both sides one obtains 



(102) 



The above expression only makes sense, if (f\{E) and 
ip-{E) are both in thus (p+{E) must be in $+. Hence 
the instructions for computing brackets while using the 
retarded unity decomposition are: On the left hand you 
have to use $+-vectors and on the right hand -vectors. 
However, since both test spaces are dense in Ti., you would 
not have to care too much, if there was not a significant 
subtlety: time symmetry is broken. To see this we inves- 
tigate the time evolution of a test vector S <&-: 



-iEt 



(103) 



Since f{E) is in i?^, there is an analytic continuation 
ip-{E — iy) to (D_ and so there is one for its propagated 
counterpart, given by 



e-y'e 



-iEt 



V-{E - iy) 



(104) 



Obviously, condition ( p8| ) is violated for and t < 
thus time evolution leads out of the test space for t < 0. 
As a consequence test functions in can only be prop- 
agated to the future. Analogically, test functions in $+ 
can only be propagated to the past. Using condition (|9q) 
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one finds in addition that the extension of the unitary 
propagator U{t) to the distribution spaces and <i>L is 
only possible for respectively t > and t <0. So retarded 
Gamov vectors \f^) in <I>^^ "propagate" (as part of the 
spectral decomposition of Hilbert vectors) only to the fu- 
ture and advanced Gamov vectors |/;7) in <i>!L propagate 
only to the past. The unitary propagator group U hence 
splits into a semigroup U+ propagating to the future (and 
applying to vectors in <&_ and <i>^^) and a second one 
propagating to the past (and applying to vectors in 
and <I>L)- This feature is called by the Brussels School 
the intrinsic irreversibility of quantum mechanics and is 
identified with a microscopic arrow of time. 



C. Time symmetry and Gamov vectors 

Is the special choice of the two RHS's necessary for the 
correct mathematical implementation of Gamov states? 
If so, one could regard the splitting of the unitary time 
evolution group into two semigroups as enforced by the 
mathematical concept. This point of view is taken by A. 
Bohm and N.L. Harshman in ^] : 

"Philosophizing alone would not be enough 
to take a semigroup instead. To arrive at the 
semigroup, we start from the empirically de- 
sirable properties of Gamov resonance states 
and let mathematics determine the path." 



The above mentioned properties are given by 



H\f+) = {En-iT/2)\f+) 
\ftit)) 



'iHt \ f+ 



(105) 
(106) 



with Efi, r > 0. Starting from these conditions, following 
the reasoning of section IVB, one is led to the two RHS's 
with broken time symmetry. However, the argumenta- 
tion is not valid if a Rigged Hilbert Space can be found, 
where Gamov vectors are included and time symmetry is 
preserved. Let us start with the space "D of C°°-functions 
with compact support. As already mentioned, this space 
is a perfect test space fulfilling the requirements given in 
section IV A. Furthermore we have V C S. Now regard 



the space Z of inversely Fourier transformed D-functions: 



Z F-^{V} 



(107) 



Since the inverse Fourier transform is a one-to-one map- 
ping on the Schwartz space S, we have 



Z dS 



(108) 



Let ip{s) be a P-function. Its inverse Fourier transform 
>f{E) G Z reads 

^(E) r dse^^^^(s) - ^ / dse^^^^(s), (109) 



where K is the compact support of <^(s). The analytic 
continuation of f{E) to the complex plane exists and is 
an entire function, given by 



tp{E + iy) = ^ I dse-y'e''^'(p{s) 
Jk 



(110) 



for aU i;, y G R (see Q, theorem 7.22). Since the in- 
verse Fourier transform is a linear continuous one-to-one 
mapping, the transformed topology T-p of P, 

Tz ■■= F-^rx,} , (111) 

transfers all topological properties from V to Z. Thus, 
{Z, Tz} is a topological vector space fulfilling all require- 
ments of a test space in L^(R). Since Z contains func- 
tions being analytic in the whole complex plane, the topo- 
logical du al contains the retarded Gamov vectors de- 
fined by (|l05| ) and (106). Furthermore, it contains also 
the advanced Gamov vectors j/^T) with eigenvalues in 
the upper half of the complex plane, as well as any other 
generalized eigenvector used in the extended spectral de- 
compositions ( |37| ) and (^J) . Because of relation ( |108[ ) all 
Z-functions are 5- functions, so Z is already a physical 
convenient test space and standard quantum mechanics 
is applicable, too. Since T) <f_ iS(IR±) <f. V the space Z 
does neither include any of the test spaces $± nor is it 
part of it, i.e. 



z (^^±(/:Z 



(112) 



but of course Z n <I>± / 0. So Z is a different test space, 
whose topological duals include any of the complex distri- 
butions needed for a comfortable treatment of decaying 
systems. 

However, there is no intrinsic irreversibility. Let the 
2?-function 'f>{s) have a compact support bounded by the 
finite interval [a, 6]. Its inverse Fourier transform ^p{E) 
can be propagated in time yielding ft{E) — e^^^*Lp{E), 
whose Fourier transform. 



dE 



-iEs -iEt 



^{E)^^{s + t) , (113) 



has a compact support bounded by the finite interval 
[a + t,b + t], thus it remains compact for all t G R and 
iptis) stays in T>, so (pt{E) stays in Z. Hence the propaga- 
tor U{t) does not lead out of the test space for any t G R. 
Gonsequently, the unitary propagator group U — {U{t)} 
can be extended to the distribution space Z^ without any 
splitting. The support of Z-functions can be restricted 
to the positive semiaxis, due to the positivity of energy, 
and the RHS 



(ZdL^C Zt) 



l]R4 



(114) 



is thus suitable for the presented concepts. Any desired 
unity decomposition including Gamov vectors, retarded 
or advanced, can equivalently be used. As far as Hilbert 
states are concerned, no unphysical behaviour is pre- 
dicted and any calculated physical quantity is numeri- 
cally identical to the one calculated by standard meth- 
ods. 
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D. Physical justification 



Apart from pure mathematics there is a physical rea- 
soning presented by members of the Brussels School who 
try to derive the special choice of the two Hardy classes 
from the assumption of causality (see e.g. Schulte and 
Twarock in j|] or Antoine, Bohm and Harshman in 
from where the citations below are taken). This is done 
the following way. 



1. Argument 

Starting point is a special assumption of causality 
which is called the "preparation— >registration arrow of 
time" : 

"Time translation of the registration appa- 
ratus relative to the preparation apparatus 
makes sense only by an amount of t > 0." 

This is translated into the language of quantum mechan- 
ics to a so-called "quantum mechanical arrow of time" : 

"A state (p^it) must be prepared before an 
observable — IV'~(0))(V'~(0)| can 

be measured in that state, i.e. = 4>'^{Q) 
must be prepared during t < 0." 

Interpreting measurement formally as a scattering exper- 
iment, the prepared state is identified with an incom- 
ing state and the measured state '0" with an outgoing 
state (see fig. H). Survival probabilities of the initial (i.e. 
prepared) state are thus calculated by transition proba- 
bilities of the incoming and outgoing state: 



P{t) 
P{t) 



'IV-")!^ fort>0, 
for t < 0. 



(115) 
(116) 



Using the energy decomposition this implies 



J dEp{E)e-'^*(j)+{E) ^0 for t > (117) 
and J dEp{E)e-''^*tP+{E)^0 for t < 0. (118) 



Using the Paley- Wiener theorem 
desired conclusion 



e H 



3) one is led to the 



(119) 



Hence, founding on the given formulation of causality the 
only correct choice of a Rigged Hilbert Space is the choice 
of the two RHS's of the Brussels School. In other words: 
Standard quantum mechanics has not yet found a micro- 
scopic arrow of time because it has been incomplete in 
its axiomatic foundation. 




prepared 
state 



measured 
state 



Past 



t = 



Future 



FIG. 9. Measurement as a scattering experiment. The pre- 
pared state is an incoming state cj)^ , the measured state is 
an outgoing state leaving the interaction region where the 
measurement takes place (grey circle). The survival ampli- 
tude of the prepared state is the transition amplitude between 
incoming and outgoing state. 



2. Critics 



The above argumentation is correct, but it cannot be 
taken as a derivation of the arrow of time. As is ob- 
vious even by the choice of words ("before", "after", 
etc.), the given "causality" assumption is already based 
on the arrow of time, hence it is logically impossible 
to derive the latter from the former. Causality itself 
does not require an arrow of time, since neither New- 
ton's law nor Schrodinger's law contradict causality, al- 
though both theories are perfectly time-symmetric. Yet, 
the "preparation— !-registration arrow of time" refers ex- 
plicitely to open subsystems, namely the preparation ap- 
paratus and the measuring apparatus, which are coupled 
together to form a composite system. Though one cannot 
generally infer from statements about the time evolution 
of open subsystems to statements about the evoution of 
closed composite systems, especially not when the sub- 
systems are involved in a measurement process. Hence 
postulating a "preparation— ^registration arrow of time" 
for all systems equals introducing a microscopic arrow 
of time by hand into quantum mechanics. Once the ar- 
row of time has been put into the mathematics, it is not 
surprising that it will reappear eventually. The situation 
can be visualized by fig. O. 



Arrow of time I 

T 



Causality | — » 2RHS's~| — 



Splitting ofU 



]- 



FIG. 10. A "causality argument" equivalent to the postu- 
lation of an arrow of time can be used to derive the special 
choice of two RHS's, which in turn leads to an arrow of time. 
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Thus, the special choice of the two RHS's must be re- 
garded as being equivalent to a postulation of a global 
microscopic arrow of time. The spht RHS's are a mathe- 
matical mirror image of the empirically motivated split- 
ting between past and future. 



E. The arrow of time as a matter of choice 

As already stated in the introduction, it is possible to 
explain an arrow of time due to time-asymmetric evolu- 
tion of macroscopic observables and due to local obser- 
vations on subsystems. But respecting a global micro- 
scopic arrow of time, there is a major difference between 
standard and Rigged Hilbert Space quantum mechanics. 
In standard quantum mechanics there is no way to im- 
plement a global microscopic arrow of time except per- 
haps by applying brute force and splitting the unitary 
evolution group axiomatically into two semigroups. In 
Rigged Hilbert Space quantum mechanics one has the 
choice of either timc-symmctric or timc-asymmctric mi- 
croscopic evolution, due to the choice of the test spaces. 
Hence, in the RHS formulation it is possible to imple- 
ment a global microscopic arrow of time in contrast to 
standard quantum mechanics. So one could argue that 
nature obviously has taken the time-asymmetric choice 
of two test spaces instead of one. But what does "ob- 
viously" mean? One might mention our everyday expe- 
riences or global time-asymmetric boundary conditions 
like the big bang or the second law of thermodynam- 
ics. Though these experiences and boundary conditions 
are all related to subsystems or far from being micro- 
scopic. The "universe" is a construction from observa- 
tions on subsystems on a very macroscopic scale. Al- 
together, those time-asymmetric boundary conditions all 
base on local macroscopic observations and hence, as al- 
ready stated, they do not give rise to a real problem 
with quantum mechanics, which is a set of microscopic 
laws. So the choice of global microscopic time-symmetry 
or time-asymmetry remains merely axiomatic and cannot 
be derived from other physical postulates. 



V. SUMMARY 

The quantum mechanical formalism, developped by 
the Brussels School, offers a comfortable treatment of de- 
caying systems, where standard quantum mechanics re- 
quires a lot more effort. The generalized spectral decom- 
positions involve Gamov states and can be interpreted in 
a physically intuitive manner. They directly yield impor- 
tant physical quantities like decay rates and energy level 
shifts. Perturbation expansions become possible in cases, 
where standard quantum mechanics fails. 

In contrast to standard quantum mechanics, where the 
implementation of a microscopic arrow of time is not pos- 
sible, the Rigged Hilbert Space formalism offers the pos- 



sibility to incorporate such an arrow by the special choice 
of two different Gelfand triplets. Due to this choice, the 
unitary time evolution group splits into two semigroups, 
each one for respectively the future and the past. How- 
ever, it is shown that the use of Gamov states does not 
require this splitting. As a counter example, one single 
Gelfand Triplet is constructed, where the formalism of 
the Brussels School can be applied and where time sym- 
metry is preserved. Furthermore, it is pointed out that 
the physical arguments in favour of the special choice of 
two Gelfand Triplets are not independent from the arrow 
of time and thus cannot be used for a derivation of the 
latter. Hence even within the framework of the Rigged 
Hilbert Space formalism, the global microscopic arrow of 
time remains a postulate. 
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